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Abstract 

A search for scalar top quarks in R-parity violating supersymmetry is performed in e + p 
collisions at HERA using the HI detector. The data, taken at y/s = 319 GeV and 301 GeV, 
correspond to an integrated luminosity of 106 pb -1 . The resonant production of scalar top 
quarks t in positron quark fusion via an R-parity violating Yukawa coupling A' is consid- 
ered with the subsequent bosonic stop decay t — > bW. The R-parity violating decay of the 
sbottom quark b — ► du e and leptonic and hadronic W decays are considered. No evidence 
for stop production is found in the search for bosonic stop decays nor in a search for the 
direct R-parity violating decay t — ► eq. Mass dependent limits on A' are obtained in the 
framework of the Minimal Supersymmetric Standard Model. Stop quarks with masses up 
to 275 GeV can be excluded at the 95% confidence level for a Yukawa coupling of electro- 
magnetic strength. 
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1 Introduction 



Deep inelastic collisions at HERA are ideally suited to the search for new particles coupling to 
an electron 1 -quark pair. Such particles include squarks (q), the scalar supersymmetric (SUSY) 
partners of quarks, in models with R-parity violation ($ p ) [1]. In most scenarios the squarks of 
the third generation, stop (t) and sbottom (b), are the lightest squarks. In the present analysis we 
focus on resonant stop quark production in eg-fusion which proceeds via an^? p coupling A'. We 
investigate SUSY scenarios in which the sbottom mass is smaller than the stop mass, M~ b < M^, 
which are complementary to previous ]/t p SUSY searches for squark production by HI [2,3]. 
This study is particularly interesting following the observation of events with isolated electrons 
or muons and missing transverse momentum [4]. The dominant Standard Model (SM) source 
for such events is the production of real W bosons. Some of these events have a hadronic final 
state with large transverse momentum and are not typical of SM W production. These striking 
events may indicate a production mechanism involving processes beyond the Standard Model, 
such as the production of a scalar top quark and its decays as proposed in [5]. 

A' ~ 

In this paper a search is presented for stop quarks which are produced resonantly, e + q — > t. 
Of particular interest is the bosonic decay t — > bW, where the sbottom decay into SM particles, 

~ A' 

b — > u e d, is also R-parity violating. This decay mode is experimentally investigated for the first 
time. The analysis includes both leptonic and hadronic W decays. A scenario is investigated, in 
which decays of the light squarks into neutralinos and charginos are kinematically not possible. 

~ A' 

In order to cover all decay modes, the $ p decay t — > e + d is also considered. The corresponding 
diagrams are shown in figure [T] At HERA, stop quarks with masses close to the kinematic limit 
of ~ 300 GeV can be produced. Such high masses are kinematically inaccessible at LEP and 
the bosonic stop decay modes considered are difficult to observe at the Tevatron. 

The analysis uses the data collected with the HI detector in positron-proton scattering in 
the years 1994-2000, where the energy of the incoming positron is E e = 27.6 GeV. The 
proton energy in 1994-1997 is E p = 820 GeV, which leads to a centre-of-mass energy of 
y/s = 301 GeV. The data correspond to an integrated luminosity of £301 = 37.9 pb~ x . In 
the years 1999 and 2000, where the proton energy is E p = 920 GeV and the centre-of- 
mass energy is s/s = 319 GeV, the data recorded correspond to an integrated luminosity of 
£319 = 67.9 pb- 1 . 



2 Phenomenology 

The most general supersymmetric theory which is gauge invariant with respect to the Standard 
Model gauge group allows Yukawa couplings between two SM fermions and a squark or a 
slepton. These couplings induce violation of R-parity, defined as R p = (_i) 3n B+ n £+ 2S \ where 
tib is the baryon number, ul is the lepton number and S is the spin of a particle. At HERA 
the Yukawa couplings between a lepton-quark pair and a squark are of particular interest [6]. 

'in the following, the term electron refers to both electrons and positrons. 
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Here the resonant production of stop quarks and the ft p decay of stop and sbottom quarks via a 
non- vanishing coupling A' 131 are investigated. Both processes are described by the Lagrangian 

£jtp ~ K31 e Lhd R + \' 131 Ve,LhdR, (1) 

where the indices L and R denote the left and right states of the fermionic fields and their 
corresponding scalar superpartners. The coupling A' 131 is a free parameter of the model with the 
subscripts 131 being the generation indices. 

In the third generation large mixings between and q R are conceivable [1]. Because of the 
structure of the squark mass matrices the stop and sbottom are the most likely candidates for the 
lightest squark states. The phases 9 q (with q = t or q = b) parameterise the mass eigenstates, 

q\ = q L cos 9 q + q R sin 9g and q 2 = -q L sin 9g + q R cos 9g, (2) 

with the convention M qi < M q2 . Since the $ p stop interaction involves only the t L component 
of the fields, the production cross sections of stop quarks scale as 

a$ x ~ Xf 31 d(x) cos 2 9i and o~ t2 ~ Xf 31 d(x) sin 2 9 h (3) 

where d(x) is the probability of finding a d quark in the proton with a momentum fraction 
x = M? /s, where 2 denotes the stop masses. The lighter state does not necessarily have 
the larger production cross section. However, in the SUSY parameter space investigated in 
this paper, Mj 2 is large enough to ensure that the resonant production of t 2 can be neglected. 
Therefore in the following the notation t will indicate the lighter t\. 

Searches for fermionic squark decays via their usual gauge couplings (into neutralinos, 
charginos or gluinos) are presented in [2]. In the present, complementary analysis the SUSY 
parameter space is chosen such that these decays are kinematically forbidden. It is moreover 
assumed that the sbottom quark b\ (denoted by b) is lighter than the lightest stop, such that the 
only possible decay modes are t — > bW with W — > //' and the HjL v decay into SM fermions, 
t — > e + d. It has been checked that the three-body decays via an off-shell W can be neglected 
compared with the ]jt p stop decay for the values of A' 131 which can currently be probed at HERA. 
Thus, only the region Mi > + Mw is investigated here, where the stop quark can decay 
into a sbottom quark and a real W. The branching ratio BR^ W for this decay mode depends 
only on the masses of the squarks involved, the coupling A' 131 and the mixing angle 9^. It 
is proportional to cos 2 9~ b . This branching ratio is shown for example values of and A' 131 as 
a function of the stop mass in figure El Under the assumption that squark gauge decays into 
fermions are kinematically suppressed, the sbottom will subsequently undergo the R\, decay 
b — > P e d and several final states can be investigated depending on the W decay mode. The 
four signatures considered in this analysis are given in table[H with the corresponding diagrams 
shown in figure[H Taking into account the LEP lower bound on the sbottom mass [7], the mass 
range chosen is 180 GeV < M~ t < 290 GeV and 100 GeV < M~ b < 210 GeV. 

The interpretation of the results is performed within a Minimal Supersymmetric Standard 
Model (MSSM) in which the masses of the neutralinos, charginos and gluinos, as well as the 
couplings between any two SUSY particles and a standard model fermion/boson, are deter- 
mined by the usual MSSM parameters: the "mass" term /x which mixes the Higgs superfields, 
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the soft SUSY-breaking mass parameters M x , M 2 and M 3 for U(l), SU(2) and SU(3) gaugi- 
nos and tan (3, the ratio of the vacuum expectation values of the two neutral scalar Higgs fields. 
The usual GUT relations between M\, M 2 and M 3 are assumed to hold [1]. These parameters 
are defined at the electroweak scale. All sfermion masses are free parameters in this model, 
as well as the squark mixings 9f and 9~ b and the soft SUSY-breaking trilinear couplings A t and 
A b . The squark mixing parameters, masses and the SUSY parameters A t , A b , tan (3 and ft are 
related by 

Ml - Ml M? - Ml 

A t = h •sin2fl t ~ + /icot/3 and A b = 61 62 • sin2fl s + //tan/3 , (4) 

with Mt and Af fe being the top and bottom masses, respectively. 



3 The HI detector 

A detailed description of the HI detector can be found in [8]. The HI detector components 
relevant to the present analysis are briefly described here. The right-handed coordinate sys- 
tem used is centered on the nominal interaction point with the positive ^-direction defined to 
be along the incident proton beam. The Liquid Argon (LAr) calorimeter is used to identify 
jets and electrons and covers the polar angle range 4° < 9 < 154° with full azimuthal ac- 
ceptance. It has an energy resolution of a(E)/E pa 12%/y/E/ GeV © 1% for electrons and 
a(E)/E pa 50%/ s/E] GeV © 2% for hadrons, as obtained in test beam measurements. The 
energy measurement is complemented by a calorimeter in the backward region [8,9]. The 
central and forward tracking detectors are used to measure charged particle trajectories, to re- 
construct the interaction vertex and to supplement the measurement of the hadronic energy. The 
central part of the detector is surrounded by a superconducting magnetic coil with a strength of 
1.15 T. The iron return yoke is the outermost part of the detector and is equipped with streamer 
tubes to form the central muon detector (4° < 9 < 171°). It is supplemented by the forward 
muon spectrometer (3° < 9 < 17°) which uses a toroidal magnetic field. The luminosity is 
determined from the rate of the Bethe-Heitler process ep — > ejry, detected in a calorimeter 
located downstream of the interaction point. The main triggers for the processes investigated 
are provided by the LAr calorimeter and their efficiencies are close to 100%. 



4 Monte Carlo event generation and simulation 

For each possible SM background source a detailed Monte Carlo simulation of the HI detector 
response is performed. All processes are generated with an integrated luminosity much higher 
than that of the data. 

To determine the contribution of neutral current (NC) deep inelastic scattering (DIS) events 
ep — > ejX, where j indicates a jet, the RAPGAP [10] event generator is used, which includes 
the Born, QCD Compton and boson gluon fusion matrix elements. Higher order QCD radia- 
tive corrections are modelled using leading logarithmic parton showers [11]. An important SM 
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background for the bosonic stop decay channels is charged current (CC) deep-inelastic scatter- 
ing, which is simulated using DJANGO [12]. QCD radiation is implemented to first order via 
matrix elements, while higher orders are modelled by parton shower cascades generated using 
the colour-dipole model, as implemented in ARIADNE [13]. In both generators QED radia- 
tive effects arising from real photon emission are simulated using HERACLES [14]. For the 
simulation of the direct and resolved photoproduction of jets, ep — > (e)jjX, the PYTHIA 6.1 
program [15] is used, which includes light and heavy quark flavours. It contains the QCD 
Compton and boson gluon fusion matrix elements and radiative QED corrections. In the above 
event generators the parton densities in the proton are taken from the CTEQ5L [16] parameteri- 
sation. The most important SM background to the leptonic W decay channels is the production 
of W bosons, calculated in leading order (LO) using EPVEC [17]. By reweighting the events 
as a function of the transverse momentum and rapidity of the W boson, next-to-leading order 
QCD corrections are accounted for [18]. The production of multi-lepton events may also con- 
tribute to the SM background for the leptonic W decay channels when one lepton is undetected 
and some fake missing transverse momentum is reconstructed. This process is generated with 
the GRAPE [19] program. 

The predictions of the RAPGAP, DJANGO and PYTHIA models are scaled by a factor of 
1.2 for cases where three jets are required. This factor accounts for deficiencies in the parton 
shower model for multi-jet production and is obtained by comparison with data [20]. 

For the SUSY signal simulation and the calculation of its cross section SUSYGEN [21] 
is used which relies on the LO amplitudes for ed — > bW given in [5]. The parton densities 
are taken from the CTEQ5L parameterisation and evaluated at the scale of the stop mass. All 
bosonic stop decay topologies are simulated for a wide range of stop and sbottom masses in a 
grid with steps of typically 20 GeV; for the $ p stop decay only the stop mass is varied. The 
events are passed through a detailed simulation of the HI detector. These simulations allow the 
signal detection efficiencies as a function of the stop (and sbottom) masses to be determined in 
the entire phase space since the mass steps are sufficiently small for linear interpolations to be 
used. The variation of the efficiencies with the coupling A' 131 when the stop mass and width are 
both large is also taken into account. 



5 Event selection and analysis 

The selection of the event topologies, as given in table [H relies on the identification of jets, 
leptons and missing transverse momentum, as detailed below. The primary interaction vertex 
has to be reconstructed within 35 cm in z of the nominal position of the vertex. Non-ep back- 
ground is rejected by searching for event topologies typical of cosmic ray and beam-induced 
background [22] and the event timing is required to be consistent with the ep bunch crossing. 

5.1 Particle identification 

The electron identification is based on the measurement of a compact and isolated electro- 
magnetic shower in the LAr calorimeter. The hadronic energy within a cone defined by R = 
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a/ (A77) 2 + (A0) 2 < 0.5 around the electron direction is required to be below 2.5% of the elec- 
tron energy. Here, 77 = — ln(tan(0/2)) is me pseudorapidity and <p denotes the azimuthal angle. 
For electron polar angles in the region 10° < 9 e < 140° a high quality track is also required 
to be associated to the electromagnetic cluster. This allows efficient rejection of photons which 
convert late in the central tracker material. 

The muon identification is based on the measurement of a track segment or an energy de- 
posit in the instrumented iron associated with a charged particle track in the inner tracking 
systems [4]. In addition, a track in the forward muon system is taken as a muon candidate. The 
muon momentum is measured from the track curvature in the solenoidal or toroidal magnetic 
field. A muon candidate should not deposit more than 8 GeV in the LAr calorimeter. The dis- 
tance between the muon candidate and the nearest track is required to be R > 0.5. For muon 
pairs a cut on the track opening angle and polar angle sum is applied to reject cosmic muons. 

Jets are reconstructed from energy deposits in the LAr calorimeter combined with well 
measured tracks using a modified inclusive k±_ algorithm [23,24] in the laboratory frame. Elec- 
tron and muon candidates are excluded from the algorithm. Only jets in the polar angle range 
7° < 6 Jet < 140° are considered to ensure that they are reliably measured. To reject elec- 
trons which are misidentified as jets, topological criteria for electron-jet separation are applied. 
About 80% of fake jets and 3% of genuine jets are rejected, as determined from simulations. 

The missing transverse momentum P± is derived from a summation over all identified 
particles and energy deposits in the event. In the channels where one or more neutrinos are 
expected, an event is only accepted if the energy and the momentum in the beam direction fulfil 
— P z ,i) < 50 GeV, where Ei is the energy and P zi is the z component of the momentum 
and the index i runs over all hadronic and electromagnetic objects and muons. This requirement 
reduces the contamination due to badly measured NC DIS events 2 where fake missing transverse 
momentum is reconstructed. 

5.2 Systematic uncertainties 

The sources of experimental and theoretical systematic uncertainties considered in this analysis 
are described in the following. They are added in quadrature. 

• The electromagnetic energy scale uncertainty is between 0.7% and 3% depending on the 
particle's impact position on the LAr calorimeter [22]. The uncertainty on the polar angle 
of electromagnetic clusters varies between 1 mrad and 3 mrad, depending on 9 [22]. The 
uncertainty on the azimuthal angle is 1 mrad. The tracking efficiency is known with a 
precision of 2% for polar angles above 37° and deteriorates to 15% in the forward region. 

• The muon P T scale uncertainty is 5%. The uncertainty on the polar angle is 3 mrad and 
on the azimuthal angle is 1 mrad. 

• The hadronic energy scale of the LAr calorimeter is known to 2%. The uncertainty on the 
jet polar angle determination is 5 mrad for 9 < 30° and 10 mrad for 9 > 30°. 

2 A NC DIS event is expected to have — Pz,i) = 2-E e = 55.2 GeV due to energy and momentum 

conservation. 
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• The uncertainty on the integrated luminosity results in an overall normalisation error of 
1.5%. 

• Depending on the SM production process different theoretical uncertainties are used. 
These amount to 15% for W production, 10% for NC DIS processes and 15% for pho- 
toproduction. For ep — > vjjjX reactions, the theoretical uncertainties are about 20%, 
which takes into account the deficiencies of the parton shower modelling of multi-jet 
production. 

• For the SUSY signal detection efficiencies, an uncertainty of 10% is assumed resulting 
mainly from the linear interpolation in the grid of simulated mass values. An additional 
theoretical systematic error, mainly due to the uncertainties on the parton densities, af- 
fects the signal cross section. This uncertainty varies between 12% at lower stop masses 
(x « 0.3) up to 50% for stop masses of 290 GeV (x « 0.8) at ^ = 319 GeV. An ad- 
ditional uncertainty of 7% on the signal cross section [2] arises from the variation of the 
scale at which the parton densities are evaluated. 

5.3 Analysis of the bosonic stop decay channels 

According to table [U the bosonic stop decay leads to three different final state topologies. If 
the W boson decays into leptons, the signature is a jet, a lepton (electron or muon) and miss- 
ing transverse momentum (jef± channel and jfi^± channel). The W decay into v T r, where 
r — > hadrons + v, is not investigated in this paper. If the W decays into hadrons the signature 
is three jets and missing transverse momentum (jjjf± channel). The selection of the final states 
analysed is described in the following sections. 

5.3.1 The channels t —> jefx and t —> jnP± 

The selection criteria for the jef± and jfxf± channels are the following. 

• A lepton must be found with P| > 10 GeV and with polar angle 5° < 6 e < 120° for the 
electron and 10° < 6^ < 120° for the muon. 

• A jet must be found with P^ et > 10 GeV within the angular range 7° < 9 Jet < 140°. 

• The total missing transverse momentum must satisfy f± > 12 GeV. 

• The difference in azimuthal angle between the lepton I and the direction of the system 
X to t, composed of all other measured particles in the event, must be A0(7 — X to t ) < 165°. 
NC background events with a mismeasured electron are rejected by this cut. 

• The azimuthal balance of the event must satisfy Vap/Vp < 0.3, where Vap jVp is defined 
as the ratio of the anti-parallel component to the parallel component of the measured 
calorimetric transverse momentum with respect to the direction of the total calorimetric 
transverse momentum [25]. This cut removes NC DIS events which generally have high 
values of Vap/Vp. 
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• In the jeP± channel, the variable y e = 1 — E' e (l — cos9 e )/(2E e ), where E' e denotes 
the energy of the scattered electron, is required to fulfil y e > 0.3. This cut reduces the 
remaining NC DIS background, since particles coming from a bosonic stop decay will 
be boosted in the forward direction, leading to a rising da/dy distribution. This contrasts 
with the steeply falling da/dy ~ y~ 2 distribution of NC DIS events. 

The stop mass cannot be reconstructed in these channels since there are two neutrinos in the 
final state. Therefore the transverse mass distributions are investigated. The transverse mass is 
defined as 

M T = yj (f j_ + P| + P} et f - {f ± +P<r + i^ et ) 2 , (5) 

where P~± , P^ and P^ et are the missing transverse momentum, lepton and jet momentum, re- 
spectively. The transverse mass distributions are shown in figures Ok and |3j). In the jeP± 
channel, 3 events are found while the expectation from the SM is 3.84 ± 0.92 events. In the 
jfiP± channel, 8 events are found while 2.69 ± 0.47 are expected. This slight excess could be 
interpreted as a scalar top decaying via t — > bW [5]. All 11 events in the jeP± and jfiP± chan- 
nels were also found in [4], where additional events were selected since there were no explicit 
jet requirements. Between 60% and 70% of the SM expectation arises from the production of 
real W bosons. The numbers of events and the SM expectations can be found in table |2j The 
stop signal efficiency amounts to typically 35%-45% for the jeP± channel and 30%-40% for 
the j[J>P± channel and depends mainly on Mi and M~ b . 



5.3.2 The channel t — > jjj f> ± 

For the jjjP± final state topology the following criteria are required. 

• Three jets must be found with P^ etl > 20 GeV, P/ et2 > 15 GeV and P/ et3 > 10 GeV, 
each with polar angle 7° < 9 Jet < 140°. 

• The total missing transverse momentum must satisfy P± > 25 GeV. 

• The selection is restricted to yh > 0.4 exploiting the different y h distributions of the stop 
signal and the SM background. Here, is calculated using y h = ^ h (E h — P Zjh )/2E e 
[26], where ^2 h (E h — P Zjh ) is calculated from the hadronic energy deposited in the de- 
tector. 

Assuming that only one neutrino is present in the event and applying the constraints P± = 
P^ and ~ Pz,i) + {E u — P z ,v) = 2E e , the neutrino four-vector can be calculated. Hence, 

the invariant mass M rec can be reconstructed in this final state topology with a mass resolution 
of about 15 GeV. In figure |3t the reconstructed mass distribution for the jjjP± channel is 
shown. A total of 5 events are found while 6.24 ± 1.74 are expected from SM processes (see 
tableEl). The SM background arises predominantly from CC DIS processes. The stop detection 
efficiency is typically 35%-50% in this final state topology. 
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5.4 Analysis of the R-parity violating stop decay channel t ed 

For stop and sbottom masses for which Mi « M~ b + Mw, the ft p decay t — > ed becomes 
dominant (see figure El). Events from this process are characterised by high Q 2 NC DIS-like 
topologies. The momentum transfer squared, obtained from the scattered electron, is defined 
by Ql = (-fr)V(l ~~ Ve)- Both the stop decay and the NC DIS final states consist of a jet and 
an electron with high transverse momenta. However, the distributions of the events in mass 
M e = ^/x e s and y e are different. Here, the Bjorken variable x e is related to the other kinematic 
quantities by Q\ = x e y e s. Stop decays via $. p lead to a resonance in the M e distribution. In 
addition, stop quarks decay isotropically in their rest frame leading to a flat da/dy distribution, 
contrasting with that of NC DIS events. 

The selection criteria for the t — > ed channel are the following. 

• The longitudinal momentum loss is limited by requiring 40GeV < ~ -^m) < 
70 GeV. 

• An electron must be found with P% > 20 GeV and with polar angle 5° < 9 e < 120°. 

• A jet must be found with P^ et > 20 GeV and with polar angle 7° < 9 Jet < 140°. 

• The total missing transverse momentum and \[P^ must fulfil P± / \[P^ < 4 V GeV, 
which takes into account the energy resolution of the LAr calorimeter. 

• Only events with Q 2 e > 2500 GeV 2 are considered. 

• The selection is restricted to y e < 0.9 to avoid the region where migration effects due to 
QED radiation in the initial state are largest. Background from photoproduction, where a 
jet is misidentified as an electron, is also suppressed by this cut. 

• In order to maximise the signal sensitivity, a mass dependent lower y e cut is applied as 
in [2], which exploits the differences in the M e and y e distributions between the SUSY 
signal and the DIS background. 

The M e spectrum for data and for the SM expectation are shown in figure 01 for all HI 
e + p data. The resolution in M e is between 5 GeV and 9 GeV, depending on the stop mass. No 
significant deviation from the SM is found. In particular, at masses above ~ 180 GeV where the 
stop signal is searched for, no significant peak is observed in the data. A total of 1100 events 
are found, while 1120 ± 131 are expected from SM processes, mainly from NC DIS events. 
The numbers of events and the SM expectations can be found in table El In the ed channel, the 
typical stop signal efficiency is about 30%-45%. 

6 Results of SUSY analysis 

6.1 Interpretation of bosonic stop decay searches 

In the jfJ,JP± channel a slight excess of events compared with the SM expectation is observed, 
confirming the previous HI analysis [4]. All other channels are in good agreement with the SM 
(see table 121). 
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Assuming the presence of a stop of mass M t - decaying bosonically, the observed event yields 
are used to determine the allowed range for a stop production cross section er f ~. The number 
of observed and expected events satisfying the relevant selection cuts, N data and N SM , are 
integrated within a mass bin (transverse mass bin) around the calculated stop mass (transverse 
mass), corresponding to the decay channel under consideration. The width of the mass bin 
is adjusted to the expected mass resolution, such that each bin contains events reconstructed 
within ±2 standard deviations of the given stop mass. A signal cross section o~ t dependent on 
the stop mass can be determined from N data and N SM in each bosonic decay channel by folding 
in the signal efficiency e, the t and W branching ratios BR^ bW ■ BR w ^jj, and taking into 
account the integrated luminosities £ 301 and £ 319 : 

r nr \ N data — N$M 1 ,~ 

a ~^ M ^ = 7~rp rd ~ — r — Tr~- (6) 

e " £>K-t-+bW ' nlt W^ff' T <J ' ^"301 + ^319 

Here, r a is the ratio of the theoretical stop production cross sections at s/s = 301 GeV and 
y/s = 319 GeV. The branching ratio for t — > bW is assumed to be BRi_^~ hW = 100%. The 
uncertainty on the cross section, Acr t -, is determined from the statistical error on the number 
of observed events and the systematic uncertainty on the SM prediction. The bands in figure 
|4|represent the allowed cross section regions for all bosonic decay channels. The band for the 
333 fi. channel is narrow due to the large branching ratio BRy/^ q q>- 

From figure|4]it can be seen that the stop interpretation of the excess seen in the jfxf± chan- 
nel is not supported by the other decay modes. For instance, the probability that the observed 
event rate in the jjjf± channel fluctuates upwards to produce at least the number of events ex- 
pected on the basis of the signal in the channel is around 1%, depending slightly on the 
stop mass. Hence, exclusion limits on the R\ v SUSY model described in Section|2]are derived. 



6.2 Exclusion limits in the MSSM 

The results from the selection channels considered in this paper are combined to derive con- 
straints in the MSSM. For a given set of parameters, the full supersymmetric mass spectrum 
and the branching ratios of all stop and sbottom decay modes are calculated. An upper limit 
Gum on the stop production cross section is calculated at the 95% confidence level (CL) using a 
modified frequentist approach based on likelihood ratios [27]. 

Each considered channel contributes via its branching ratio, the signal efficiencies and the 
number of observed and expected events within sliding mass bins (transverse mass bins). Al- 
though the selection criteria for the various channels are not explicitly exclusive, it was checked 
that double counting of events is negligible. The given set of model parameters is excluded if it 
predicts a cross section which is larger than cr iim . 

In order to investigate systematically the dependence of the sensitivity on the MSSM para- 
meters, a scan of the SUSY parameter space is performed. The SUSY parameter space is 
selected such that the combined branching ratio is 

BR tot = BR- t ^ ed + BR^- bW ■ BR- b ^ Ued > 85%. (7) 
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The parameter M 2 is set to 1000 GeV and fx is restricted to the range 400 GeV < ji< 1000 GeV, 
which ensures that the gaugino masses are large. The mixing angles 0i and 6~ b are allowed 
to vary between 0.6 rad and 1.2 rad. For given values of tan/3, A t and A b , the parameters 
, Mi , 9f, 6i and \i are scanned. Here, A t and Ab are only needed to determine the masses of 
the heavier stop and sbottom, according to equation ©, they are set to A t = At, = —100 GeV. 

For each point in the 5-dimensional parameter space an upper bound on the coupling A' 131 
is obtained. The resulting limits are given for two cases: (i) every point of the scanned SUSY 
parameter space is excluded, (ii) at least one point in the scanned SUSY parameter space is 
excluded. The resulting limits obtained for tan j3 = 10 are shown in figure |5k and \5]p in 
the (Mf, Mi) plane for A' 131 = 0.1 and A' 131 = y/4ira em = 0.3. At A' 131 = 0.1 stop masses 
Mi ^250 GeV can be excluded, while masses Mi £ 275 GeV are excluded at a Yukawa coupling 
of electromagnetic strength, i.e. A' 131 = 0.3. The resulting limits projected on the (M t ~, A' 131 ) 
plane for M~ h = 100 GeV are shown in figure|5k. For Mi = 200 GeV, couplings A' 131 ^0.03 are 
ruled out and for Mi = 275 GeV the allowed domain is A' 131 ^0.3. The limits do not signifi- 
cantly depend on tan (3 or on M 2 , provided that M 2 > 400 GeV, which has been checked by 
repeating the analysis with tan (3 = 2 or M 2 = 400 GeV. 



7 Conclusions 

A search is performed for scalar top quarks resonantly produced in e + p collisions at HERA 
in R-parity violating SUSY models. Final state topologies resulting from R-parity conserving 
bosonic stop decays or R-parity violating direct decays are considered. In the jfif± channel, 
a slight excess of events compared with the SM expectation is observed. Nevertheless, no 
evidence for stop production is found, since this excess is not supported by the other three 
channels analysed in the present paper. 

For the first time, direct constraints on stop quarks decaying bosonically are derived. Includ- 
ing the direct $ p stop decay, mass dependent limits on the coupling A' 131 are obtained within 
a Minimal Supersymmetric Standard Model. In a large part of the model parameter space, the 
existence of stop quarks coupling to an e + d pair with masses up to 275 GeV is excluded at the 
95% CL for a strength of the Yukawa coupling of A' 131 = y/Aira. em = 0.3. 



Acknowledgements 

We are grateful to the HERA machine group whose outstanding efforts have made this ex- 
periment possible. We thank the engineers and technicians for their work in constructing and 
maintaining the HI detector, our funding agencies for financial support, the DESY technical 
staff for continual assistance and the DESY directorate for support and for the hospitality which 
they extend to the non DESY members of the collaboration. 



13 



References 

[1] For reviews see: 

H. P. Nilles, Phys. Rept. 110 (1984) 1; 

H. E. Haber and G. L. Kane, Phys. Rept. 117 (1985) 75. 

[2] A. Aktas et al. [HI Collaboration], "Search for Squark Production in R-Parity Violating 
Supersymmetry at HERA", submitted to Eur. Phys. J. C, DESY-04-025, hep-ex/0403027 

[3] C. Adloff et al. [HI Collaboration], Eur. Phys. J. C 20 (2001) 639 [ |hep-ex/01 02030] . 

[4] V. Andreev et al. [HI Collaboration], Phys. Lett. B 561 (2003) 241 |hep-ex/0301030| |. 

[5] T. Kon, T. Matsushita and T. Kobayashi, Mod. Phys. Lett. A 12 (1997) 3143 
| |hep-ph/9707355| |. 

[6] J. Butterworth and H. Dreiner, Nucl. Phys. B 397 (1993) 3 | |hep-ph/9211204 |. 

[7] LEP SUSY working group, http:/lepsusy.web.cern.ch/lepsusy/. 

[8] I. Abt et al. [HI Collaboration], Nucl. Instrum. Meth. A 386 (1997) 310 and 348. 

[9] R. D. Appuhn et al. [HI SPACAL Group Collaboration], Nucl. Instrum. Meth. A 386 
(1997) 397. 

[10] H. Jung, Comput. Phys. Commun. 86 (1995) 147. 
[11] T. Sjostrand, Phys. Lett. B 157 (1985) 321. 

[12] G. A. Schuler and H. Spiesberger, Proceedings, Physics at HERA (1992) vol. 3, 1419. 
[13] L. Lonnblad, Comput. Phys. Commun. 71 (1992) 15. 

[14] A. Kwiatkowski, H. Spiesberger and H. J. Mohring, Comput. Phys. Commun. 69 (1992) 
155. 

[15] T. Sjostrand, P. Eden, C. Friberg, L. Lonnblad, G. Miu, S. Mrenna and E. Norrbin, Comput. 
Phys. Commun. 135 (2001) 238 |hep-ph/00100"T7| . 

[16] H.-L. Lai etal, Eur. Phys. J. C 12 (2000) 375 |hep-ex/9903282|. 

[17] U. Baur, J. A. Vermaseren and D. Zeppenfeld, Nucl. Phys. B 375 (1992) 3. 

[18] K. P. Diener, C. Schwanenberger and M. Spira, Eur. Phys. J. C 25 (2002) 405 
| |hep-ph/0203269| |; 

P. Nason, R. Ruckl and M. Spira, J. Phys. G 25 (1999) [hep-ph/9902296]; 
M. Spira, PES Y-99-060, |hep-ph/9905469l 

[19] T. Abe, Comput. Phys. Commun. 136 (2001) 126 |hep-ph/0012 029|. 

[20] S. Caron, "Jets in photoproduction at HERA", RWTH Aachen, Germany (2002), DESY- 
THESIS-2002-035. 

14 



[21] S. Katsanevas and P. Morawitz, Comput. Phys. Commun. 112 (1998) 227 
| |hep-ph/97 11417| ; 

N. Ghodbane, S. Katsanevas, P. Morawitz and E. Perez, 'SUSYGEN 3', hep-ph/9909499 



[22] C. Adloff etal. [HI Collaboration], Eur. Phys. J. C 30 (2003) 1 [hep-ex/0304003|. 

[23] S. D. Ellis and D. E. Soper, Phys. Rev. D 48 (1993) 3160 |hep-ph/9305266|. 

[24] S. Catani, Y. L. Dokshitzer, M. H. Seymour and B. R. Webber, Nucl. Phys. B 406 (1993) 
187. 



[25] C. Adloff etal. [HI Collaboration], Eur. Phys. J. C 13 (2000) 609 [hep-ex/9908059]. 

[26] A. Blondel and F. Jacquet, Proceedings of the Study of an ep Facility for Europe, ed. 
U. Amaldi, DESY report 79-48 (1979) 391. 

[27] T. Junk, Nucl. Instrum. Meth. A 434 (1999) 435 [ |hep-ex/9902006l . 



15 



Channel 


Decay process 


Signature 




t -> bW 






'—t dv e 




jef± 


W -> ez/ e 


jet + e + f ± 




— > tv t — > ez/z^i/ 








jet + /i + f ± 




— ► TV T — > jJLVW 




jjjF± 


W -> gg' 


3 jets + f j. 


ed 


i ^ ed 


jet + high P T e 



Table 1 : Analysed stop decay channels in $ p SUSY. The $ p processes are indicated by the 
coupling A', and f± denotes the missing transverse momentum. 



Channel 


V« = 301 GeV 


y/s = 319 GeV 




combined 




data 


SM expectation 


data 


SM expectation 


data 


SM expectation 


jef± 


1 


1.16 ±0.28 
(W: 0.75 ±0.12) 


2 


2.68 ±0.64 
(W: 1.80 ± 0.29) 


3 


3.84 ±0.92 
(W: 2.55 ±0.41) 


mfi. 


4 


0.84 ±0.14 
(W: 0.57 ±0.09) 


4 


1.85 ±0.33 
(W: 1.36 ± 0.22) 


8 


2.69 ±0.47 
(W: 1.93 ±0.31) 


333 f± 


1 


1.91 ±0.54 


4 


4.33 ± 1.21 


5 


6.24 ± 1.74 


ed 


366 


384 ± 45 


734 


736 ± 86 


1100 


1120 ± 131 



Table 2: Total number of selected events in the various stop decay channels for the HI e + p data 
at y/s = 301 GeV, </s = 319 GeV and the combined data set. For the jef± and channels 
the SM expectations arising from W production are given in brackets. 
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a) 

Figure 1 : Lowest order s channel diagram for 
bosonic decay of the stop and b) the ft p decay 




b) 

\ p stop production at HERA followed by a) the 
if the stop. 




Figure 2: Examples of the stop branching ratios as a function of the stop mass for 
Mi = 100 GeV and A' 131 = 0.1, when the fermionic decay modes of the stop via the usual 
gauge couplings are kinematically suppressed. The solid lines show the branching ratios for 
9~ h = 0.6 and the dashed lines for 9^ = 1.2. The sum of the branching ratios is slightly less than 
one since hadronic r decays following W — > vr are not considered. 
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Figure 3: Mass spectra for the HI e + p data: a) transverse mass of the jeJP± channel; b) trans- 
verse mass of the jfif± channel; c) reconstructed mass of the jjj f± channel; d) invariant mass 
distribution of the ed channel. The data are compared with the SM expectations with the system- 
atic uncertainties shown as the shaded band. The expected signal from a t with mass 260 GeV 
is also shown with arbitrary normalisation. 




260 280 

M T [GeV] 

Figure 4: Bands representing the allowed stop cross section regions oi ± Acr f - as a function of 
the stop mass as obtained from the analysis of each bosonic stop decay channel. 
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Figure 5: Exclusion limits at the 95% CL in the (Mj, Mg) plane for a) A' 131 = 0.1 and b) 
A' 131 = 0.3. c) Exclusion limits at the 95% CL on the ]jt p coupling A' 131 as a function of the stop 
mass for = 100 GeV. The limits are derived from a scan of the MSSM parameter space as 
indicated in the legend. The two full curves indicate the regions excluded in all (dark) or part 
(light) of the parameter space investigated. 
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